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ABSTRACT 


The effect of shear on the rate of po 1 yme r i zat i on of styrene 
in presence of activated carbon has been studied . It was found 
that the rate of po 1 yme r 1 zat i on decreases in presence of shear 
and is enhanced in presence of activated carbon - Based on the 
chemistry, a kinetic model has been proposed in which carbon is 
participating in the po lyme r 1 za t i on and is found to describe the 
experimental results well. Ule have also carried out simulation 
of free radical polymer ization . Analytical expressions for 
isothermal po lyme r i zat i on have been developed for monomer and 
polymer c on cen t rat i ons and for moments of radical and dead 
polymer distributions in terms of a dimensionless variable Z. 

The analytical solution has later been adopted for non i s o t he rma 1 
polymer izat ion in presence of gel effect . 
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CHAPTER i 


Introduction 

The chain growth polymerization has been studied in detail 
and has been shown to consist of initiation, propagation and 
termination steps. These can be schematically written as 
Initiations 


21 


I+M U P * 


Propagation: 


Termination: 


n 


-► P 


n+1 


P'+ I ► |y| 
n n 


p*+ p* — ^ hr 
n m n+m 


k t . 

P*+ P" M + M 

n m n m 


( 1 ) 

As can be seen above, the initiation is a sequential 
reaction of decomposition of initiator molecules to give primary 
radicals and subsequent combination of primary radicals with 
monomer molecules. In reality all primary radicals do not 
produce propagating chains and they remain in 'cages* consisting 
of surrounding molecules before they diffuse apart. However 
before this diffusion takes place, the primary radicals may react 
with one another to form an inactive species. In order to 
account for this recombination, a factor called initiator 
efficiency, f is introduced. There are several factors which 



affect the Initiator efficiency and these have been discussed in 

14-18 

the literature in detail 

In the propagation reaction, the monomer adds on to the 
polymer radical and the propagation rate constant k^ has been 
assumed to be independent of chain length. Polymer radicals may 
be converted into polymer molecules by primary and mutual 
(combination and disproportionation) termination. The expression 
for rate of polymerization is obtained by making use of the 
standard kinetic assumption of a steady state approximation of 
transient species, namely, the chain radicals. Thus we obtain. 


r = k 
P P 


r f k d a 2 ] V 

|_ <k tc + k td ) . 


1 

2 


CM3 


( 2 ) 


This equation is valid only for low conversions and it will 

be shown later that the Quasi Steady State Approximation <QSSA> 

breaks down due to the gel and glass effects. 

- n 

A parameter y indepe/dent of monomer and initiator 
concentration may be defined as 


y 


T72 


CM-ICIgU 


(3) 


The rate of polymerization at low conversions. in the 
absence of shear for several monomers have been determined 
experimentally and excellent compilations are available^. 
Studies have been carried out establishing the effect of mixing 
upon the rate of polymerization in which simple stirring devices 
were used. However, no correlation exists which relates the rate 
of polymerization with shear rates. In this work we have studied 



the effect of shear on the rate of polymerization of styrene 

using azo iso butyronit rite (AIBN) as the initiator- We have 

also studied the effect of activated carbon of various particle 

size on the rate of polymerization. 

It is well established that after certain conversion- the 

'gel effect’ appears as a sharp increase in the rate of reaction, 

usually accompanied by an increase in molecular weight and 

viscosity. The gel effect occurs after about 10% conversion 

depending on the temperature and the amount of initiator used. 

Numerous theoretical and experimental studies of the gel effect 

have been reported and it is generally agreed that it is due to 

the reduction in the termination rate constants, k t and 

. After a conversion of 70-80% is reached, the reaction 

mass becomes highly viscous and the chain propagation becomes 

diffusion controlled. As a result of reduction in k , the rate 

P 

of polymerization goes down as the reaction temperature 

approaches the glass transition temperature of the polymerizing 

29—33 

mixture. This behaviour is known as the glass effect 

Most of the studies of chain growth polymerization in 
reactors assume isothermal reactor operation. In commercial 
practice however, reactors cannot be maintained in isothermal 
conditions. Practical problems of heat transfer and thermal 
runaway invariably occur Because of their simplicity and 

low costs, tubular reactors are commercially employed and in the 

literature they have been theoretically as well as experimentally 

♦ ^ 35-42 
investigated 

In earlier studies of free radical polymerization. it was 


common to assume that the quasi steady state approximation holds. 



4 


Under this approximation, the rate of formation and rate of 

depletion of the polymer radicals are equal. Therefore, 

k^CPD 2 = 2fk ,CI„D <4) 

t s d 2 

Here CP3 s represents the total concentration of polymer 

radi cals with QSSA valid and is the initiator concentration. 

However, when gel effect sets in, k^ has been shown to decrease 

in value, which will result in increased CPU under QSSA, this 

s 

way leading to higher rates of polymerization. Since all 

polymerization reactions are exothermic in nature, the 

temperature of the polymer mass would go up for the same cooling 

rate which in turn would give higher rate of initiation. Mence 

under the influence of gel effect, the upward thermal drift 

always occurs and the QSSA breaks down. 

If the QSSA breaks down, the mole balance relation for 

polymer radicals CPU would be governed by a nonlinear first order 

differential equation instead of the simple relation given in eqn 

(4). The concentrat ion CPI in the reaction mass is normally a 

very small number compared to the monomer or initiator 

concentration. As an example, for styrene polymerizing at 60°C 

with benzoyl peroxide initiator. CPU^ calculated from eqn (4) 

—8 9 

is of the order 10 moles per liter. As a consequence in the 
numerical computation of CPU through its governing differential 
equation, the time increment At must be chosen very small. This 
would mean that the differential equations governing the 
performance of the tubular reactor carrying out radical 

polymerization are numerically extremely 'stiff 1 and the 
numerical solution is computationally time consuming. 

In this work, we have written down the transient mole 



balance relation for all apeclea in the reaction maaa. Aaauaing 
various rate constants for the mechanism of isothermal free 
radical polymerization given in eqn (1) as time independent, we 
develop the complete analytical solution. We have also 
integrated the energy balance equation making use of the above 
assumption. After every time step At we estimate the 
temperature, and rate constants are evaluated at the new 
temperature which are made use of in calculating CMD, CP3 and 
other quantities in the next time step. This removes the 
stiffness of the differential equation and results are obtained 
conveniently in negligibly small computation time. 



CHAPTER 2 


Theoretical Davelopments 

It is possible to make mole balance for various soecies in 
the reaction mass and the batch reactor results are given in 

Table 1. We define concentration of polymer radicals, CPU as 

oo 

CPU = E CP 1 (5) 

n=1 n 


With the help 

expression for the 
dCPZJ 


~Sl 


of eqn (1.2) of Table 
time variation of CPU 
= 2fk d CI 2 3 - k t CPU 2 


1, 

as 


one can 
below 


derive an 


( 6 ) 


When the gel and thermal effects are present in radical 
polymerization, the rate constants kj and k^ are dependant on 
temperature and monomer conversion and therefore cannot be 
integrated. We will show that eqn (6) has a solution for 
constant k^ and k^. and then develop an analytical solution of 
isothermal polymerization in absence of gel effect. Subsequently 
we will show that these results can be naturally adopted in the 


presence of the gel and thermal effects. 

Sol ut 1 ori of" Isothermal Radical Polymerization* 

The mole balance for the initiator given in 

Table 1 can be integrated for time invariant fk^ 

-fk t 

CI„D = CI„: e 
c. So 

In the above eqn, CI„D is the concentration of i 

do 

Let us now transform CPU in eqn (6) 


eqn (1.1) of 
as 

(7) 

nitiator at t=0. 


CPU = 


1 dy/dt 

-£7 ~V~ 


(S) 


and upon substituting eqn (7) for CI^U in eqn (6) one obtains, 


2 f k , t 

= 2fk .k t CI„3, e } 

dt 2 d t 2 o 


( 9 ) 



Further, we define x related to time of polymerization 

-f k .t 


aa 


This gives 


* “ z,k d l t [ V. • 


dx 


■art — fk d- x 
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dy _ dy dx 

tTt ~ clt = ~ fk d x 


d2 y - d dy dx 

^ > m 


dy 

3x 


=<fk d > 


2 X -Sr? + ( f k j ) 2 X 2 


lx • w *d 
In terms of these, eqn (9) becomes 

( f k . ) 2 x „ . Wli 
d .2 ' ' ^d 

we further transform x by 


dx* 


(b) 

( c) 

<d> 


2 2 d^v o j 

— h + (fk > 2 « d y _ 

ttk d > x ^ xy = O 


where 


Z = 2Vbx 

b = kty 


(a) 

(b) 


in 


terms of which eqn (12) reduces to 
r Z rf2 y . , dy 2 

Z ^ z sz- ~ (Z ~ o > y=o 

For this v has the following solution 

y = c 1 i o (z.> + c 2 k o ( 2 ) 

Where I Q (Z) and K q (Z) are the modified zeroth order 
functions and C ^ and C 2 are the constants of integration, 
the help of eqn (8), one can get CPI as 

f k . Z 


CPI = - 




1 dy 
y m. 


fk d z - C 1 I 1<Z> 

2k t ^o (Z)+ t: i JL o (Z) 


( 10 ) 


( 11 ) 

( 12 ) 


(13) 

(14) 

( 15) 
Bessel 
With 

( 16 ) 



(17) 


fk d Z K n (Z) - C I^Z) 

= 2FT K <Z> + C 1777“ 

to o 

where C = C^/C^r which is to be determined by the initial conditions 

Let us assume that at t=0, the concentration of Polymer 

radicals is CP3 . In the Z plane, time t=0 corresponds to Z 
o u 

given by 


Z = 2 
o 


A 


2k t CI 2 3 o 
f k 


(18) 


and C in eqn ( |7 ) is given by 


C = 


2k CP] 

K.(Z ) - -FT— ° K (Z ) 

1 o ffc .Z o o 

d o 


2k CPI 

W + TFT7 ” W 


d O 


K.(Z ) 

1 o 

TTTZ y 

1 o 
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„ K ( Z ) 
* o o 

o T7TT 7 
1 o 


1 + CPD 


7 I (Z ) 

* O Q 

o TTTZ ) 

1 o 


( 19) 


where 


* CP:, o 

CPI1 o = rF37 - ^ 


. / 


2fk 

d 2 o 


( 20 ) 


Above CPH s is the concentration of polymer radicals assuming 
the steady state approximation and is given by eqn (4) 

The variation in the monomer concentration, CM1 , is governed 
by eqn (1.3) of Table 1. We substitute eqn (13) in it to obtain 


dCMU 

~ <rr~ 


dCMT 

~a t~ 


dz 

~sz 


f k .1 

a 


dcr-o 

dz 


k CMDCPD 
P 


or 


f k ,Z 
d 


dCMT 

Sz~ 


fk Z 

k p CM Vr 


K , ( Z ) 
1 


C I^Z) 


TTT7T 

o 


TTTTZT 

o 



d[M] 

dz 


( 21 ) 


9 


k r k cz) - c i^z) " 

CM3— j— — K (Z) + C I (Z) 

t L o ° J 

We substitute 

u = K (Z) + C I (Z) 
o o 

which on differentiating with respect to Z gives 
ijH =- K ^ ( Z ) + C I t (Z) 

Comparison of equations (19) and (21) gives 

— dCM3 k p du 

"TFO IT 

which on integration leads to 


CM3 = CM3 

o 


K (Z) f C I (Z) 1 k p /k t 
o o 

K (Z ) + C I (Z ) 
o o o o J 


( 22 ) 


(23) 


(24) 


(25) 


Above CM3 is the monomer concent rat ion at t=0 (or Z=Z ) and Z 
o o o 

is defined in eqn (IS). 

The magnitude of Z as defined in eqn (13) is very large. 

For example, for methyl methacrylate polymerizing at 60°C with 

AIBN CI 0 3 = 0.0258 moles/litl initiator , the rate constants k,, 

2 o d 

—3 —1 5 

k and k* are 0.475x10 min , 0.4117x10 lit/mol min and 

P t 

lO 

0.20383x10 lit/mol min respectively. Taking initiator 

efficiency f to be 0.58, we get at time t=0, as defined in 

eqn (13) and x q defined by eqn (10) to be 1.3175x10^ and 

4 

2.8756x10 respectively. Therefore Z at time t=0 is given by eqn 

(18) and is equal to 0.1231x10 7 . Hence the asymptotic expansion 

of Bessel functions valid for large arguments (Z > 5) may be 

.8 

used . 

Neglecting terms involving Z in the denominator, we get 


= K (Z) 
o 



K . ( Z ) 
1 


e 


Z 



I.(Z> - I <Z) 
1 o 


/2nl 


( 26 ) 


Making use of the above approximations for Bessel 


functions, C defined in eqn (19) gives 


C = y n e 


where y 


1-CP3 


1+CP3 


(27) 


The eqn. (17) for CPU may be written as 


fk .Z 
d 


K ^ ( Z ) 


I^Z) 


kTT 

o 


After substituting for C and Bessel functions, we get an 
expression for the time variation of CP3 as 


CP3= 


fk .Z 1 - ye 
d 


-2 ( Z -Z) 
o 


(28) 


1 + ye 


-2 ( Z -Z) 
o 


The quantities CI o (Z) and CI o (Z q ) in eqn (25) for CM3 are given 


-Z -(Z -Z) 

o o 

Cl (Z) = yne e 

o r 


-/2*Z 


Cl (Z )^ yne 
o o 


(29) 


These mav be taken as zero, since both Z and Z are of the 

o 

order 10^*. Hence the eqn (25) for CM3 simplifies to 


CM3 * CM3 e 
o 


-m ( Z -Z) 
o 


( 30 ) 


k 

_ P 

m - -j— - 

K t 


where 



By making use of the above expressions for [M] and [P], we 
can obtain analytical expressions for various moments of radical 


and dead polymer distributions. 

Radical and Dead Polymer Moments: 

The differential equation for first moment G for radical 

_ <*J 1 

distribution [ Cp n 3 J can be derived from eqn.1.2 of Table 1 

as follows 

dQ 1 

ST" =_k t CP;iQ 1 +2fk d CI 2 :i+k p i:,vi;i{2i:p 1 ^ + 3CP 2 D+4CP 3 D+ > 

-k p CMDCCP 1 3+2CP 2 2+ > 

= "»< t .CPDQ 1 +2fk d CI 2 D+k p CM3{(1CP 1 D+2CP 2 3+3CP 3 D+ ) + 

(CP 1 3+CP 2 3+CP 3 3+ )} 

-kpCMTf CP 1 3+2CP 2 3+3CP 3 3 > 

dQ 1 

'SfT = 2fk d CI 2 3-k t CP3Q 1+ k p CM3CP3 (31) 

The above equation is a linear differential equation of 
first order if we assume rate constants as time invariant. We 
get the integrating factor (I.F.) as follows 

I-F = exp tk t J CP3dt> (32) 

Observing that 

d-CCMU/CMIl > 

TLMJ/CmA = - k p cp:] dt (33) 

eqn (32) reduces to 

1 - F = exp { - -i ln < rm >} 

= fCM3 \~ Um 

\CM3j" (34) 

The integration of the differential equation (eqn 31) 


Q -| as 


gives 
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sOSt}"' “ J <2fk d C V + V M cp:,> {rraJ 

Substituting Cl^l from eqn (13) in terms of Z as 


dt 


2 

Z c fk 


rt j = 

LJ S J 8k 

and carrying out the integration to finally obtain 

v s[Sy- 2 [-y 1/n, + ^_<z +1) 


where 


- CM3 

C 1 = 1-1/m and c 2 = Q 


CM3 f k . 
o d 


10 T^T7aT 21T 


(Z +1) 
o 


(35) 


(38) 


In order to derive the differential equations governing 

zeroth, first and second moments of dead polymer distribution, we 

make use of the mole balance for M from Table 1. We have 

n 


dCM 3 
n 

dt 


k + n-1 
t c 


J 

m— 1 


co 

LP n? LP „ n? + K A ZP 1 T CP ' 
m n— m td • n A . n 

n-1 


(37) 


Further we know that. 


oo n-1 

XT* Lf 
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y n k y 

CP 3CP 3 = Y CP 3 \ 

n-E m-1 

m n-m n /. 

n-1 m-1 


The zeroth moment of dead polvmer distribution X. — 

o 


be derived by adding eqn (37) for all n 


dX k 

sz ~ =-ir- LP1 + k td CP:, ‘ 


where CP3 


and 


= y cp 3 

n-1 n 

eqn (38) is made use of to obtain eqn (39) 
Making use of eqn (33) 

-1 


CP3 dt = 
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o 


~TT TTFD7TFnTT 
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we get. 
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(39) 
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I 


CP3 dz 


( 40 ) 


Substituting for CP3 from eqn. (28), we get 

fk 

p ~ L t 1+we o 


2(Z -Z) 
o 


-X = JdX = -(0.5 k + k + .) ~ fZ. 1 ~y > - 

o oo J o tc td k_ 2kT J “ -2(Z -Z) 


dz (41) 


The rigorous evaluation of the integral leads to an infinite 
series and is given in Appendix 1. However, if Z above is 
assumed to be approximately equal to Z q , the integration in 
eqn (41) is considerably simplified without affecting results 
very much. On doing so eqn (41) leads to 

fk jZ 


X - X 
o o 


— f dX = -(0.5 k„ 
o J o t 

This can be easily integrated to 
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The differential equation for first moment 


fx = ' 

n CM3 1 

l 1 l 

L n J 


n=1 


of dead polymer distribution is similarly derived using eqns (37) 
and (38) as follows: 
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On eliminating [P] with the help of eqn. 1.3 of Table 1, 


_ 'H 

mr =--ir p 


dtCMH/OI ]> 
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Substituting for and using 


diCMD/CMD > 
o 


m dZ 


We get. 


x i x io 


CM3/CM3 


.md2 


where X 1Q is the value of X jJ at time t=0 and C 1 and C R are as 
defined in eqn (36) 

On integrating, we obtain. 


X 1 X 10 


- ~ 8 (i-©y ' 1 ^ 


(45) 


The expression for second moment of Polymer radical 

2 

distribution E n CP^U) is dervied making use of eqn (1.2) of 

Table 1, to obtain 
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This is again a linear differential equation of first order 
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Q r™2_i m = rpo 'i 

2 U M V J imp 


Jr 

m 


2fk CI D+ 2k CM ft:M3 -> 
n 2 P o [CM3J 


Q.+b CM 3CP3 
1 p o 


dt+ constant 


Substituting for Cl 1 ^nri n 

y ‘-■ l 2 J and Q i' we get 

1 

*-y 

m 


- f£M3 ") m r z fk 

McMIiJ “ J 2fk d ~8k^ 


z fk d <z n~ z) O 

d e ° 2-dZ , 

^fk + | 2 V« 




+c 


1 

m fk 


fCM3 T ‘ tk d 

amp + 2k^ 


(Z+1 ) 


■ 

1 


_ 2-dZ. r 

ey 5 -! 

rCM3 ^ 


c fCM3 1 


+c * (47) 


We integrate the above e vnr a . • 

expression assuming Z=Z . j^ e 

rigorous solution has also been Wc - . . ° 

derived and is presented i n 

appendi x X and we find t k a t ♦« 

we find that the r, suUs derived in the following 

gives results within 0.01%. 


fk .(Z+1) 

0 o = d 

2 2k~ 


4k CM3 C 
P Q 

< .Z (2m — 

d o 


fCMJ - 

. 2 4k CM3 C_ 

_ _ P o 2 | 

fCM3 'J 


> f k Z m 

d o 

[CM3 J 

*• or 


i+i 

m 


- iSSitay- 

..... C* I. .......... ... . >t i-c 

Therefore we get 


C* = Q 


20 


fk d (z 0 +i) 

2JT 


4k CM3 C 
P o 1 


4k p CM3 o c 2 2k^CM3 


+ Tk d Z o (km - i) + T~k d 2° m + 


2k D CM:3 n CM3 

+ T- P 0 ^ O 

F tV m-1) /m 


(49) 


The differential equation for second moment X g =£ n 2 CM □ of 



dead polymer distribution la derived u m , . 

derived using Mole balance relation 

(eqn 37) and eqn (38) as follows: 


‘td cp]Q 2 + k -¥ y> 2 y c p n ^ cp 3 

— * i i n 


k td cp:,Q z + -p 


k td CP * IQ 2 + 


X CP n‘ ] X (n>2 + ° 2 + 2nm) CP 3 
L —‘ n 

(CP3Q 2 + CP3Q 2 + 2q 2 , 


(k tc + k td ) + K_ Q 


2 *tc U 1 


(50) 


We observe that 


K «•? 


dt = 


J[ C i [cm:J + C 1 (my * (z z + 2 zai) 

* ** t&s* 1 - + '1 0$[ sy(z + l] 

f k — 

+ Cp K) ( z ♦ 0 q « 91 , 


Integrating the above equation assuming Z = Z ( 
and taking it out of the integral, we get 


in the denominator 


"to °1 dt - 


.fk _ 2, 


fk I [ z "ft i 1 - [Sy ^ + f 1 - (rSy *j 

f k 2 

+ fcEf) [ 2 ( Z 0- ZZ ) + E (V*) «"(£)] 

+ [i - [sy "*)♦ £ ^ (i - [ay] 


2C 1 C 2 
(m+1 ) Z 


* ^ ft (' - (fsy) • (■ - (ss-f) 

* ^(’ - (sy‘] 1 


( 52 ) 



17 


Also, 




CP3 Q dt = — — ! 
2 k 


r 2 

,.fCM3 3 

f“a_l 

d lcM3^J 


Substituting Q_ from eqn (48) anrf 

2 - nd carrying out the integration 

assuming Z = Z q in the denominator we get 


I 


k t CPDQ 2 dt = ~ 
P 


mfk 

i 

2k, 


2 ( Z »- z 2 ) + M] 


2k p CM3 o C 1 <• 
,K d^o'‘ ! “ ” 11 


’-(syl 


4k CM3 C 
P o 2 

f k Z (m + 1 ) 
do 


2k CM3 


1+- 

„ fCM3 3 m I 2k CM3 f ^ 

rN ]-y^(ysy) 


m ° 


1 , 


+ - c * Msy 

v j 


(53) 


where C is as defined in eqn (49) 
The r ef ore 

,2 


2k 

X 2 ~ X 20 = ~fk 


K [ f 1 - (ray] - (™y"j 

■fk H 

+ k) [ -a K-y * 

+ ^rr; [’ - (rSy "] + ir -rf f 1 - (rSyl 


+ ^ ( y f 1 - (rayj * -k f 1 ■ (ray’] 


o'J 

1 
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CHAPTER 3 


Computational Techniques 

We have derived the analytical solution for isothermal 
radical polymerization in batch reactors. In order to check the 
consistency of these results, we have computed the various 
reactor parameters using 4th otder Runge Kutta numerical 
technique and checked our results against the analytical 
so lut ion * 

The flow chart of the fourf order Runge Kutta technique is 
given in Figure 1 and we have varied At in Table 4 and seen its 
effect upon the numerical solution. For stable solution, we find 
that the time increment At must be about 0.0001 min. The flow 
chart of computation for our analytical solution is given in 

Fig. 2. It must be observed that our analytical solution are 

0 

given in terms of Z variable. For isothermal reactor there is a 
straight forward exponential relation. However for nonisothermal 
cases, there is an ambiguity in terms of properly defining Z. 
This arises because in general, the rate constants vary with time 
of polymerization. To overcome this difficulty we recall that in 
the Runge Kutta technique, the total time of polymerisation is 
divided into small incremental time. The concentration 

parameters, CM3, A q , A^, A^, CP3, and for the next 

incremental time are calculated using their values at a given 

\ 

time. This fact has been shown in fig.^. In such computations, 
between the incremental time the reaction parameters are assumed 
to be time invariant. 

For the cases when the temperature of polymerisation as 
well as gelation, we similarly divide the total time of 



Initial values 
M , P , I2 ,Qi ) Q2 
A 0 1 Ai > X 2 1 T 
Parameters 

k P , , f , Cp , A H r ,T W , dia 


t = 0 


Use RK method to 

Solve balance eqns (Table 2) 

d [I 2 1 /dt eqn. (2-1) • dCtydt eqn-(2-5) 

d[P]/dt eqn. (2-2) , dA-|/ dt eqn.(2-6) 
d[M]/dt eqn. (2-3 ) ; dQ 2 /dt eqn.(2-7) 
dAo /dt eqn. (2-4 ) dT/dt eqn (2-9) 

dX 2 /dt eqn. (2-8) > 

1 


Update variables estimate 
rates constants at new temp 


t = t + h 


if / 0 


Out 

ig.1 Runge Kutta Scheme tor tree radical polymerisation 









Fig. 2 Analytical scheme tor free radical polymerisation 
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polymerization into email incremental times. Between these time 
increments, we assume that the rate constants are time invariant. 


If we assume this all the results for isothermal reactors are 
valid, between the incremental time and only the energy balance 
equation remains to be integrated, which we do now. 

The modelling of non isothe rmal polymerization can be done by 


inclusion of the following energy balance, 

pC = ( —AH ) k CM3CP3 - (T-T ) 

p dt r p D w 


(55) 


Here, (— AH^) is the heat of the reaction, U is the heat 
transfer coefficient and T y is the surrounding temperature. 

The above eqn. may be rearranged using eqns. (1-3) and (30) 
and may be written as follows; 


where 


dy . -mx 

+ fty = «e 

y = T-T 


w 


X = z -z 

o 

ft = 


8U 


i3p"c Tk~T 

p d o 


(a) 

(b) 

(c) 


(56) 


(—AH )CM3 .m 
r o 


(d) 


ft* 


(57) 


The integrating factor of eqn. (56) is e and using the 

initial condition T=T at x=0, we obtain the following solution. 

o 


ote 


— m(Z -Z) 
o 


T-T = 

w 


(ft—m ) 


_iO / -r 71 ~ft ( Z — Z ) 

-ft (I -Z) ote o 

+ (T -T ) e . ° - 7-5 r — (58) 

o w ift-m) 


This result has been built in the algorithm in figure 2. 
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CHAPTER 4, 
EXPERIMENTAL 


(a) Chemicals: 

9 

Analytical grade styrene (Burgoyne Burbidge & Co.) is washed 
with 1 .ON IMaOH two to three times to remove the inhibitor 
hydroquinone . In order to remove traces of NaOH left behind, 
styrene is washed thoroughly with distilled water. The washed 
styrene is stored over anhydrous calcium chloride atleast for two 

days in a defreezer, before it is used for our experiments. 

/ 

The initiator used is 2-2'azo isob utyronitri le and was 
obtained from SAS chemicals & Co., Bombay. For further 
purification, it was dissolved in methanol and recrystallised 
before use. The activated carbon obtained from S.D. fine 
chemicals, was subjected to dry grinding to obtain various size 
ranges .Parti cles of size range -250 mesh and -350 mesh were used. 
Di 1 aioMior i 

A sketch of special dilatometer which combines the features 
of an ordinary dilatometer and rheometer is shown in figure 3. 
The instrument consists of a glass cylinder C (I.D. 48.6 mm) 

having a precision bore capillary A (4.45 mm) exactly at the 
centre of the cylinder. A stainless steel rod of uniform 
thickness (o.d. 3.9 mm) passes through the capillary A and 
connects the rotor R (o. d. 42.5 mm) to a suitable motor. The 
rotor rests on a glass pivot which facilitates an easy rotation 
of the rotor. The precision bore capillary B (I.D. 2.0 mm) 
attached on the side of the cylinder C is used for filling the 
monomer solution in the annular space. The calibration of the 

4 

dilatometer has already been done . 




Mercury 



For exact aeasur amenta , the volume of the dilatometer has to 
be determined at the temperature at which polymerization is to be 
carried out (i.e. at 60°C). However, since the coefficient of 
linear expansion of glass is more than 2 orders of magnitude 

(49) 

smaller than the coefficient of cubic expansion of the 
liquid to be used, the change of the volume'' of the dilatometer 
with temperature can generally be neglected. 

Procadurei 

The monomer is taken in a flask and nitrogen is bubbled 
through it for about 3 minutes. Vacuum is applied to remove any 
residual nitrogen. This way we ensure that no dissolved oxygen 
is present in the reaction mass. The dilatometer is filled with 
mercury upto height H so that the reaction mixture fills only in 
the annular space between the rotor, and the glass vessel. The 
reaction mixture consisting of monomer, activated carbon and 
initiator with known concentration is charged to the dilatometer 
upto reference mark as shown in fig. 3. The dilatometer is 
subsequently placed in a constant temperature bath which is 
maintained at 60°C. Immediately, the levels of the liquid in the 
capillaries A and B begin to rise due to the expansion of the 
liquid but subsequently it falls because of polymerization. The 
change in liquid level with time in both the capillaries A and B 
were measured using a cathetometer . The change in volume of the 
reaction mixture with time is calculated. The rotor R can be 
rotated at a constant speed by a suitable motor. The speed of 
rotation is measured by a stroboscope and is maintained at a 
constant value by carefully monitoring it. After polymerization, 
the dilatometer is emptied in the following manner. The 



26 


dl latoaet er la Inclined carefully over a small beaker so that the 
mercury can run off, which is then washed with benzene and 
acetone and finally dried. 

The rate of polymerization is calculated by standard 

di latometr i c techniques . The shear rate is calculated using the 

rheometric equation, 
w R, 

m 

r = 


R 

r g“ r r 


(59) 


where, is the angular velocity of rotor, R_ and R,_ are radii of 

K vj 

rotor and cylinder respectively. 

The dilatometric equation is 


y 


p Ps* p 


m 



1 

T72 

cmdci 2 t 


(60) 


where, Pp s and p m are densities of polymer in solution and 
monomer respectively. M^is molecular weight of monomer, V is 
the volume used, CMH and CI^D are the concentrat ion of monomer 
and initiator respectively. 

The experiment is conducted for different concentration of 
activated carbon and also for different size range of activated 
carbon . 
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CHAPTER S 

Results and Discussion 

O m 

The earlier works of Kumar and Shukla have shown that the 

changes observed in the rate of polymerization, r of 

P 

sytrene— benzoyl peroxide system with shear, y is purely due to 
variation of initiator efficiency. 

The polymerization of styrene has been carried out in a 
special dilatometer as previously described. The change in 
height in the capillaries has been measured and converted into 
the change in volume of the reaction mass. Regressional analysis 
of experimental runs has been carried out and slopes of AV 
versus time plot are used to calculate the rate of polymerization 
Tp. The parameter y as defined in eqn (3) has been calculated 
and results have been summarized in Tables 7,8,9 and 10. 

We have carried out polymerization of styrene under 
following conditions: 

(1) Polymerization of styrene with AIBN initiator in the absence 

of activated carbon. *••• •" ' '' / 

(2) Polymerizat ion of styrene with AIBN initiator in presence of 

* » ”, . -,. w , -v 

activated carbon (100 mg/200 ml styrene). ^ ' r 

(3) Polymerization of ,sty rene— AIBN system in presence of 
different particle sizes of activated carbon. 

In order to study the effect of shear on rate of 
polymerization, parameter y has been defined, which is 
independent of monomer and initiator concentration. This has 

been plotted in fig. 4. From theoretical considerations, we 
should have y independent of shear rate. Contrary to this, we 
find that with increase in shear rate y first goes down and then 



With Carbon 
Without Carbon 
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Fig. 4 Variation of y with shear rate for Styrene AIBN system 
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increases to about 95% of its initial value. Ue have conducted 
experiments to show that it is the initiator efficiency which is 
affected by shear. 

The most general kinetic model for the 5 cage effect has been 
4 

discussed by Shukla and is given below: 

^COO-OQC^ - > C2<£C00 . □ 

C2^C00.D ► C^C0O?!.+C0 2 D 

Ca^COO.D+M *■ <^C0 2 . + <ji>COOM. 

LP&COQ.l ► 2<£C00. 

<£C00. + M ► ^COOM. 

0COO. — - ■ * 4 > .+ COg. 


4 > . + M *• . 

4>.+ <jt>COO . ► 4>Z0O4> 

2.<p * 4 , ~4 / 


The square bracket indicates the presence of cage and <p 
represents phenyl groups formed in the decomposition of benzoyl 
peroxide initiator. Instead of considering this detailed model, 
we simplify and consider the following three idealization in the 
initiation step. 

<1) The initiator decomposes reversibly in the cage to give 
initiator radicals and N 2 and the radicals add on to monomer 
to give polymer radicals. 

I 2 — - » 21 + IM 
I. + M ► P 

(2) The initiator decomposes reversibly to give primary radicals 
from which N 2 is eliminated subsequently and both primary 
and secondary radicals add on to monomer to generate polymer 


radi cals. 
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V N 2 


I + M 


I.+ M 
1 


(3) The same sequence of steps as above is assumed, the only 
difference being that the initiator decomposition is assumed 
to be irreversible. 


V N 2 


I + M 


I i+ M > p 

In appendix 2, these kinetic schemes have been considered 

and analytical expressions for the rate of polymerization have 

been derived. If we assume that the reactivities of monomer with 

* 

primary radicals are comparable to k and those of k , A and k,. 

P ci • dP 

« 

are comparable to k^ and k^j, it can then be shown that there is 
negligible effect of Ng concentration upon f. This suggests that 
models on cage effect alone are unable to explain the phenomenon 
experimentally observed. The fall in the rate of polymerization 
(and a subsequent rise) suggests that the Ng gas dissolved in the 
reaction mass must be acting as a catalyst on the formation of 
primary radicals as follows: 



Ng catalyst 


21 


The rate of polymerization can then be derived as 


r = k 
P P 


| fk d t=l V [:i ?_ =r j 


1/2 


CM3 


On the application of shear the concentration of Ng falls 
due to enhanced mass transfer. As a result the rate of 

polymerization falls as shown by above equation. However. for 
large shear, the cages present initially are destroyed and the 
initiator efficiency recovers as has been observed experimentally. 

We have carried out experiments in presence of activated 
carbon and we find that the rate of polymerization is increased 
in its presence. We notice that the activated carbon, has a 
distribution of particle size and we decided to separate them 
into —350 and —250 mesh screens. Subsequently we activated them 
by putting them in oven for 4 hours. In figure 5, we find that 
the smaller the particle size the larger is the rate of reaction 
and hence y Q value. We also examined the polymerization of 
styrene in presence of carbon as well as the shear rate. When 
shear is applied the rate of polymerization falls down. In 
presence of carbon, when shear is applied a similar behaviour is 
observed except that rate of polymerization is higher in presence 
of carbon at a given shear rate. 

9 

There have been several studies involving activated carbon 
in the past and these studies focussed on fixing styrene polymer 
chain to carbon particles. The affixing of styrene chain on 
carbon particles can occur only when carbon particle was 


participating in the process of polymerization. Based upon the 
this we propose that when a given carbon particle is bombarded 
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Fig. 5 Variation of y o with activated Carbon Concentration for 
Styrene AIBN-system 



33 


with an initiator radical it forms some kind of an activated 

■Jf. 

species CC 3 . We pr o pose that this activated species can react 
with another car bon to form an inert molecule or can react with a 
monomer to give a grafted po 1 yme r molecule. In presence of 
car bon , the mechanism of polymerization can then be written as 
Initiation : 


I + C 


* T 

C + C — 


Inert (C-C) (c) 


* A 

C + M — - 


1+ M 


Propagation 


P+ M 


Termination 


P + P 
n m 


(61 > 


The derivation of rate of polymerization in presence of 
carbon has been given in appendix 3, and the following expression 
has been derived. 


r — r 
p po 


1 * 

nrr + K 2 


(62) 


where K* and K* are constants dependent on temperature, 
monomer and initiator concentration as follows 


k 2 + k^CMH 

~v ~~D — c 

K 1 K 4 K d 
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* 

K 2 = 


TTF71T7 
14 d 


CI 2 3CM3 


(b) 


(63) 


For the two particle size activated carbon, we have carried 

this exercise in Fig. 6. We find that we can well represent the 

data by straight-line according to eqn.(62) and the constants 

and Kg, are different for different particle sizes. 

We have considered the batch reactor simulation of radical 

polymerization. The basic differential equations governing this 

are given in Table 2. The numerical integration of radical 

polymerization is difficult for the following reason. It may be 

recalled that the molar concentration of * intermediate polymer 

—8 

radicals is of the order of 10 . As opposed to this the monomer 

concentration is of the order of 10 while the reactor temperature 

2 

is of the order of 10 . Due to wide variation in the ranges of 
various reactor parameters, the differential equations of Table 2 
are numerically 'stiff* in nature and require selection of very 
small time increment At. In order to remove this numerical 
stiffness, we have attempted in this thesis, to determine an 
analytical solution of radical polymerization in batch reactors. 

We have assumed the reactor to be isothermal and trasformed 
the mole balance of polymer radicals into Z domain and have 
successfully obtained the analytical result. Since the value of 
Z is very large it was possible to develop an approximate ion in 
eqn.(28). The testing of this relation with those of Runge-Kutta 
numerical reuslt are given in Table 5 and is seen that both 
results are almost identical for the entire range. 

For nonisothermal case, we divide the total time of 




Fig. 6 Plot of 
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polymerization Into email Incremental time. Between thla 
incremental time we assume the reactor performance as isothermal 
and use the analytical solution of the reactor developed so far. 
The integration of the energy equation has been done in eqn . (58) 
and the algorithm has been described in figure 2. With this 
algorithm, the results have been computed and summarized in Table 
6. In the same Table the numerical results obtained from 
Runge-Kutta technique are also given. We find on comparison that 
both results are close to within IX in the entire time range. 
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CHAPTER 6 
CONCLUSIONS 

In this work, we have studied the effect of shear on the rate 

of polymerization of styrene. It was observed that r^ decreases 

on application of shear and then partly recovers. We have 

considered several models of initiator efficiency and on careful 

examination we find that none of these . can explain the 

experimental observation in the entire range.lt is then suggested 

that N^evolved in the decomposition of initiator has a catalytic 

effect on the addition of monomeric units. The application of shear 

helps removal of from the reaction mass with a consequent fall 

in r . At high shear rate there is a likelihood that the evolved 
P 

Ng, is unable to disengage from the reaction mass. As a consequence 
the rate of polymerization goes up due to higher concentration of 
Ng in the reaction mass. From the past studies, it is known that in 
presence of activated carbon , the latter participates in the 
initiation step of polymerization. We have consequently developed a 
kinetic model in which carbon participates in the polymerization. 
The model suggests that the specific surface area of carbon is 
important , which would in turn mean that small particle size of 
carbon would give higher activity. We separated commercially 
available activated carbon into -250 and —350 mesh particles and 
we find that the latter exhibits higher activity , this way 
confirming the proposed kinetic model. 

Simulation of free radical polymerization has also been 
carried out. The nonlinear differential equation for CPU has been 
reduced to a form which has its solution in terms of Bessel 
function, by a suitable transf ormat ion . Analytical expressions 
have been derived for various moments for the case of isothermal 
polymer izat ion and this has been adopted for nonisothermal 
polymerization in presence of gel effect. 
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Appendix I 

Rigorous Solution of Dead Polymer and Radical Moments 

We have assumed that Z = Z Q in eqns.(41) (47) and (52) in 

the derivation of equations for X , Q and X„. The rigorous 

o 2 2 

solution -For the above distributions, is given here . 

We find that in the derivation of the equations for dead 
polymer and polymer radical distributions, the integral. 



arises naturally. 

Making use of the t ransf ormat ion Z -Z == u, 

o 

we can write 


(A 1.1) 


r -au 

= -JlT^7 

J o 


du 


r oo 

1 ST 1 r " -au . 

- z ) l ' du 

°L £=o z o J J 


Since 


I u " 


au , au 

e du = e 


n 

V (-1 ) n ! u 

/ , ..r+1 

L — ‘ ( n— r ) ! a 

r=0 


we get the integral as 

n 


1 = e 


—a 


n n! (Z -Z) n_r 


(Z -1) sr 1 V" o 

0 X, 7 n+i X i 

z_ (n-r)! 


r+1 


r=0 o r=0 


The zeroth moment X is therefore given by 

o 


(A 1.2) 


(A 1.3) 


(A 1.4) 


K - — 1 


(0.5k* +k #J ) mf k 
c td d 


2k k t 
P t 


I 


. —2 ( Z -Z) 

1-ys e o 

-2<Z -Z) 
1+y e o 


dZ + Const 


We can write 
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• , — 2( Z -Z) 

7 e ° H7 

— 2( Z -Z) 
l+V' e o 


consider the integral. 


w 


. -2(Z -Z> 

° „ <1<Z -z> 
e' 2(Z o- Z> 


r , -2( Z -Z) 

+|(Z -Z) ^~ v e d(Z -Z) 

J 0 e~ 2 <Z o- Z> 0 


J‘V 


. — 2( Z -Z) 

Z) 17- * ° d(z -z, 

i4v <r zlz <r z > 


= x dx 

J 1+V e 

This can be integrated by parts by assuming 


1 — y e 

1+*t> e 


I t 1 i^i !r zx 1 dX 

V = X + In + v 1 * e 

Jvdu = J^X + In + v* e ^ 


However, we know that 


[ In (a + 

J X 


(A 1.5) 


cannot be expressed as a finite combination of elementary 
functions^. Therefore we express the term involving logarithmic 
function as a series 


In (1 + Y) 


=1 


n 


where Y = e 
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Therefore 


f «2 r r , 2 — 4X 3 -6X - 

| Viu -Vft* -54-*^ )dx 


00 


_ X ^ V" (-1) V 

2 ^ Z? 

n=1 


n _n 0 ~2nX 


(A 1.6) 


Hence we get 
— 2X 


00 


1 


X 1 »' dX = *f + X ln(1+v e - ZX ) - £ r: 2n - 


(A1.7) 


1+y/ e 


n=1 


It is also seen that 


: -f 


— 2(Z — Z ) 

Z I - ? ° d(Z -Z) = -Z V 

—2(Z_— Z) o o 


j . (. \ i. 

1+V" e o 

■ - z »( 
writing 


X + ln(1 + y e 


m f k 


(A 1.8) 


A = (0.5k. + k . ) p- 

tc td 2 k _ k + 
P t 


we have 


[• 


<0 


_ A| .,,2 , 2 ,., , —2 ( Z - Z ) lA V ^<“ 1 > V 

X =A | — ( Z — Z )+Z In ( 1+y e o )+) — 

n=1 


n n — 2n ( Z — Z )' 
e o 


o 1 2 o 


2n 


+C 1 (A1.9) 


where is the constant of integration to be evaluated us»ing the 

initial condition, X = X at t=0 (or Z — Z ) 

000 o 


We get. 


as 


C = X 
1 00 


Hence we finally obtain 


x -x 

o 00 


(A1. 10) 
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The equation for second moment of polymer radical 
distribution (from Eqn (46)) is 


—.-1 * 

Q 2lTOrJ = J [cmTj + 2k p CM:I Q 1 + k p CM;I CPD dt 

This can be split into three integrals, as 


Q, 


[tSy 


1 

m 


k CMUCPldt 
P 

(A 1.11) 


The first integral 1^ gives 
fk d 

I 1 = 2k~ (Z + 1} (A 1.12) 

t 


The second integral Ig is 



interms of (Z q -Z) and making use of eqn 

(A1.4) and dividing by integrating factor, we obtain for the second 
integral 


Exressing (gg-) 



*2 - 


4 V M Vi -2.(2 

fk J 6 


CO 


n 


°- Z) 


n! (Z - Z) 
o 


n— r 


— i 7 

n=0 o 


r ^ Q (n— r ) ! (Em— 1 ) 


r+1 


4k p CM:l o C 2 — (m+1 ) (Z -Z) 

— fir — e ° 


CO 


n 




n! (Z -Z) 
o 


n -r 


n=0 Z o 


r ^Q (n— r) ! in 


r+1 


2k CM3 

E 


M3 r 
o 

t 


— (m) ( Z -Z) 
e o 

(m-1) 


2k CM3 
P o 


00 n 7 .n-r 

Cm) ( Z -Z)V" 1 V " o ~ Z) 

o ) — ) -r* CA1.13) 

r n Z_j /- - v r+1 


L 


n=0 Z o 


Similarly the third Integral I 3 is 


*3 " 


CM=, o f CM3 "1 

1- 1 I™ 3 J 

m 


(A 1.14) 


For the sake of conveneience, let 
2 to 


I ‘>'-1.5*1 


" n ! C Z -Z) n r 
o 


HSo 2 o ' ^ 0 <n - ^,!!<r ' , ' 1 


(A 1.15) 


Then 


fk 

°2 = 2iT 


(Z+1) - 


4k CM3 C. _ 
p o 1 — 2m( Z 


fk 


4k CM3 C„ , „ w , 

p o 2 -Cm+1 ) (Z -Z) 

^ e o 

tK d 

2k CM3 


°- z T 
I <»> - 


( 2m- 1 ) 
2k CM3 


E - ° e—<V Z >V < m -. 

m (A 1.16) 


1 

m 


L 


Qj, = Q„„ at Z = Z 
2 20 o 

„ fk . 2k CM3 CM3 

• C 2 = Q 20 “ 2E7 <Z o + 1) + rTTnPTT + “ — T 

t t 1— — 

m 


(A 1.17) 


The differential equation governing second moment of dead 


polymer distribution (eqn 50) is 

dX„ 


dt 


t 

t c 1 


k^CP3Qg + k 



4 


Integrating eqn (51) using eqn (A 1.4) and eqn (30) we 


obtain on substituting the limits 

j-»i . r >i 2 , 


- _ 2k^ 

k .c 8 i dt = -Firi 


iSf] {s [ Z ^- z2 ) + 2[z o - Z ] - ln(z Q -z] ■ 


Further 


_ .2 -2m (Z — Z ) V , , ,.2 — H( Z -Z)V\^ 

^ e o ) <2m) C £ e o ) (2) 


CPDQ 2 dt 


D J 


d) 1( 1-e 

t J m 

-m(Z -Z) 
□ 


<Z -Z) , 
o )-' 

' °E 

d r cm i 

(CM3 1 

Um 3 q J 

d UmJ 

°2 ( 
CM3 m ( 
CM3 q 

CM3 1 A7 
CM3 J dZ 

, dZ 





(A 1.18) 


Integrating eqn CA 1.16) using eqn (A 1.4) we get, 
fk T. . „ ~ , .1 3CM3 . . -i 


: ZP2Q Z dt “ 5k 


S[i - 


3CM3 , . — m( Z — Z ) . 

o (1-e o ) 

(m-1 ) 


2k CM3 
P 

k. 


Ml 

o — m( Z -Z)V? . 

e o y Cm- 1 , i 


4k CM3 C 

L P 0 


w 

-^72m-1,J 


4k p CM3 o C 2 -<m + 1)(Z -Z)V7 -CZ-Z). 

e o )(m,m+1 >+C 2 < 1-e o ) 

CA 1.19) 


Where C* 


is as defined in eqn (A 1.17) and 



4 £ 


oo 




fl ! 


n-r ( n— r ) f ( Z -Z) n - r - p 
o 


r r+1 


n"^0 Z o ' 1^0 (r, " r)!X ^0 


Z \ ii” r / ; i 

n (n-r- 


p>! Y p+1 


(A1.20) 


Therefore, 


2k 


X 2 X 20 


t c 


fk 


, 2 — 2 ( Z -Z)r 

(2m) - Cg e o y 


2 

( 2 ) 


+C 


-m(Z 


[(ft {? [M + 2 ( Z =- Z ) + KV Z )} 

i 

2 -2m (Z -Z) V" 

i e ° L 

f fk cH 1( 1-e“ m<Z o~ Z> )_ f fk cT| 

lL 2 k t J m 

C 2 ft<--'V Z >,-C 2 (^e-<V Z > £>] 

1 orMi — 

4 (M-M 


o Z) ^T<m) 


fk 


3CM3 — m( Z -Z). 

o(1-e o ) 

(m— 1 ) 


4k CM3 C, _ 

, p o 1 -2m(Z -Z)\7-, ^ \ 

+ ! -?t e o )(2m-1,2m) 


fk 


Lu 


4k CM3 C , „ w _ 

p o 2 -(m+1)(Z — Z )\ . . „ . 

£n e o )(m,m+1) 

tK d 




jr—% | — i m ( Z Z ) 3 

2k CM! e o 


P o 


S"" 1 ’ 


. «# , , — (Z — Z ) . 

m) + Cg(1-e o ) 

(A 1.21) 



APPENDIX 2 


Various kinetic models 

The bulk polymerization of styrene is carried out in 
presence of 2,2*-azoisobutyronitrile (AIBN) having the following 
structure 


N = N 


2 r 2 * azo i sobuty roni t r i le 

This is stable at room temperature, but decompose thermally 
above 40°C giving substituted alkyl radials and nitrogen. The 
radicals remain in a ’cage* before they diffuse apart, and during 
this time a recombination of radicals may take place. In order 
to account for this, the initiator efficiency *f has been 
introduced. 

We have modelled the cage effect in three different ways to 
explain the effect of shear upon initiator efficiency *f* and 
finally show that the evolved has a catalytic effect in order 
to explain the decrease in rate of polymerization of styrene on 
application of shear. 

MODEL 1 

The following scheme is proposed: 

Initiation: 


\ 


21 + N, 


I + M 



50 




(A 2.5) 


However, we have from eqn (2) 


1 



Hence, on comparing the above two equations, we get 

M 013 

2k d CN 2 :3 + k i CM:3 


f 


(A 2.6) 
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MODEL 2 

Initiation 


d \ 


21 


s * 


T dl \ 

1 = . . ^ + N 

^ k' 
dl 


I + M ► P 

k 2 

I^+ m £* P 

Propagation 

k 

P + M ^ P 

Termination 


"RART 

■J3 




P + P 


■+ M 


d 


Balance on CI3 gives 
dCI3 


(A2.7) 


dt" " k d C1 2 3 " 2 k d CI3 


k dl C13 


+ k d.1 CN 2 :] “l 3 


k 1 CI3 CM3 (A 2.8) 


Balance on gives. 


dL 13 

- = k d1 CI3+k d1 CN £ 3 CI.,3 - k g CI,,3 CM3 (A2.9) 


- d.3- 


k a id3 


k ai CN 2 3 + k 2 CM3 


(A 2.10) 


. CI3 


k CI„3-2k ,CI3-k,,CI3+k .CN„3 — 
2 C3 dl 12./ 


k ^CI3 


k dl CN £ 3 + k 2 CM3 


- k 1 CI3 CM3 = 0 

k d CI z 2 


(A 2.11) 


2k d + k cH + k 1 ~ ~ 


kj . k,.CN_3 

di 2 


k ai CN 2 3 + k £ CM3 


Now 



Let 

Then 


52 


dCP3 

~~dX~ 


CP3 


k 1 HID CM3 + k 2 CI.,3 DID! - k t CPU (A 2.12) 

2 


= |k 1 CI3 CM3 + k 2 CI^ Cmj . 


k CM3<k .CI3CM3+k_CI <| 3CM3> 
k p CM3CP3 = -E 2 T7 g — =— 2 


(A2. 13) 


k d<\ k A<l ^P 3 

= 2k + k + k CM3 ; Q - - P 1 £■ (A2.14) 

k 4l CN 2 3 + k 2 CM3 


k^ CI3 CM3 = k 1 CM3 . k^ CI 2 3 


k 2 CI 1 3 CM3 = k 2 CM3 . k 1 CI3 


kj 1 CN £ 3+k 2 CM3 


k 2 k dl CM3 k^ CI 2 3 


k^ 1 CN 2 3+k 2 CM3 X 


(A 2.15) 


. r =k CM3 
P P 


f k d c Vf 

' ' 

k.CM3 k_k . „CM3 

1 _ + 2 41 

l k t J ] 

X X ^kj 1 CNg3 + k 2 CM3J 


(A2.1&) 


Comparing the above equation with eqn (2) we get 


k . CI3 

f - — — + 


k E k dl Cm 


( k dl 


(A 2.17) 


CC0 2 3 + kg CM3 1 
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Model 3 


Initiation: 


-* HI 


X=^ I i + N 2 

^ k' 
k d1 


I + M 


-*■ P 


1,+ M - 
1 

Pr opaga t ion: 


-*• P 


P + M 


-+ P 


Ter mi nation: 


P + P 


-► M 


(AH. 18) 


From CI<j 3 balance. 


CI^ = — 


k ai CD 


k^CNgU + k £ CM3 


From CIU balance. 


Kj cim 
cm = -S-Y - 


(A H.19) 


where 


Y = K dl + k 1 CM3 - —j 


K di K di CN 2 3 


(A H. HO ) 


k_j 1 CN 2 3 + k 2 CM3 


From CP3 balance 

k t CP3 2 = k 1 CI3 CM3 + k £ CI^ CM3 


k t CP 3 


k„k .CI 0 3CM 3 
Id H 


+ k 2 CM3 


K ai K i CI 2 3 


Y [ k dl CN 2 3 + k 2 CMIi ^ 


CP3 = 


^C! g y kl 


CM3 


K d1 K 2 


’Odl 


CN-3 +k„CM3 
H H 




J 


(AH. HI) 


r = k CM3 CP3 
P P 



54 


= k CM! 
P 


k 1 CM3 

^ Y 


K d 1 K 2 CCO 


V Y K 1 CN 2 ^k 2 CM 3 ] 

Comparing the above equation with eqn (2) we get 


l 


1 

2 


<A2.22> 


f 


k 1 CMD 

Y 


+ K dLl K 2 CI13 

Y ( k i CN 2 :,+k H i:M:i ] 


(A 2.23) 
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APPENDIX 3 

Model to account for the effect of activated carbon 
on rate of polymerization 

It has been found that the rate of polymerization goes up in 

the presence of activated carbon and the increase appears to 

2 3 

depend on the specific surface area (ft /ft ) of carbon particles 
used. The following model is put forward to explain the effect 
of activated carbon concentration on the rate of polymerization; 
Initiation: 


2 

I + C 


-► 21 
k 


1— * C* 


S k„ 


4 

C + C 


-♦ Inert (C - C) 


* 

C + M 


-4 P. 


I + M - 
Propagation 


-+ P, 


P + M 


-4 P 


Termination : 


P + P 


n m 


-► M 


(A 3.1) 


Above C is an activated species formed by the reaction of 

primary radicals with carbon. Assuming a steady state hypothesis 

(S.S. hypothesis) for C , we get 

k^CCDCII - k_ CC*3 - k„ CC*3CC3 - k.CC*3CM3 = 0 
12 3 4 


or 


„ k. CC3 HID 

* 1 

CC 1 = k 2 + k 3 CCI + k 4 CMD 


(A3. 2) 



5 < 


Balance on CID gives (assuming S.S. hypothesis for 
intermediate species) 

CIgD - k i HID CCD + k 2 CC*D - k 5 CID CMD = 0 

Further if we assume that the interaction of initiator 
radical and carbon particles is small in comparison to 
decomposition of initiator and formation of CP^D we get 

k . cim 

C13 “ T^ncra- (A 3 - 3) 

The rate of polymerization is then given by 

r = k. CC*D CMD + r (A 3.4) 

p 4 po 

Where r is the rate of polymerization in absence of 
activated carbon. 

From eqn (A3. 2) we get 


k. k . CCD C1D CMD 
4 1 

r p r po k 2 + k 3 CCD + k 4 CMD 

1 k 2 + k 4 CMD k 3 

T—=-J— = kV k: CCD CID CMD + kV k A TTT CM3 

p po 1 4 .14 


On substituting for CID from eqn (A 3.3) we get 




6 


CMD 


CIgD CCD 


k 1 k 4 k d 


CIgDCMD 


(A 3.6) 


If only CCD is varied, with initiator and 


monomer 


concentration remaining almost the same, we can write 
1 


K 


r - r 
P Po 


1 •* 
CCl + K 2 


(A 3.7) 


k„ + k. CMD 
2 4 

k A k . ki 
1 4 d 


* 

where 


(a) 
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K 


* 

2 


k„ k- •< J 

14 4 


ci 2 3cm 


(b) (A 3.8) 


A plot 



1 

versus 


is found to be 


linear and 


is 


included under results and discussion 



58 


Table 1 

Mole and energy balance relation for various species in 
batch (or tubular) reactors carrying out Radical Polymerization 


(1) Initiator, Ig 
dCIgH 

-HI— = " ^d^Z 3 


( 1 . 1 ) 


(2) Growing Polymer radi cals , CP^H , of chain length n (n=1 ,2,3. . . . ) 

t» 

CPU =f CP ] 
n=1 n 

dCP 3 

— = kj-CinCMD- k p CP 1 3CM3 - k^P^CPH 
dCP 3 

-3T- = k p CM:, [ cp „-i HP n 3 ] ' k t CP n 3CP3 M * 2) 

3. Monomer, M 
dCM3 

-3TE - XpCPDCMD (1.3) 


4. Inactive Polymer, CM^D, of chain length n (n=2, 3, . . . . ) 


dCM 3 
n 

~dl 



n— 1 

r CP 3CP 3 
. m n-m 
m=1 


00 

k t .CP 3 E CP 3 
td n n t 1 n 


(1.4) 


5. Energy balance 


p C - (—AH >k CM3CP3 - ^ (T-T ) 
p dt — r p D w 


(1.5) 
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Table 2 


! Differential Equations for fr 


ee radical polymerization 


dCIgH 


= - ^ d ci 2 3 


( 2 . 1 ) 


dCP3 


2 fk d CI 2 :3 ~ k t c:p ^ 


( 2 . 2 ) 


dCM] 

Si = ~ kpCPUCMI] 


(2.3) 


-31 - (0 - 5k tc + k td >cp:|2 


(2.4) 


2fk d C1 Z 3 ~ k t cp;i <*-, - 


kpCMICPI 


(2.5) 


~Zl = <k tc + k td )CP:iQ i 


(2.5) 


2f k d CI 2 D + k p CMl(2Q 1 + CPD ) - |< t CP3Q ( 2 .7) 


= (k tc + k td >CP:iQ 2 + k t c Q ? 


( 2 . 8 ) 


(-AH )k CMUCPU 
r P 


4U 


(2.9) 
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22,43 


Table 3 s Rate constants and other parameters for PMMA 

V 2 


A 

= 

0.168 - 8.21 x 10 6 

<T- 

B 

s 

0.03 


C 

P 

= 

0.4 cal/gm°C 


E d 

= 

30.66 Kcal/mole 


E P 

= 

4.35 Kcal/mole 


E t 

= 

0.701 Kcal/mole 


E 

Y 

%p 

= 

4.09 Kcal/mole 

28.2 Kcal/mole 


E &t 

= 

34.4 Kcal/mole 


f 

= 

0.58 


AH 

r 

= 

13.9 Kcal/mol 

AL. _ -4 

6.32 x 10 min 


■h 

= 

5.88 x 10 9 lit/mol 

min 

k° 

P 


2.95 x 10 7 lit/mol 

min 

k° 

Y 

= 

3.956 x 10" 4 


T gP 

= 

387 K 


p 

= 

5.4814 x 10 ^min 



= 

1.1553 x 10 227 mol-m 


m 


CM3 


0.9665 - 0.0011 (T-273) gm/cm" 

, _ . 3 

1.2 gm/cm 

0.0258 gmoles/lit 

8.9334 gmoles/lit 
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Table 4 


Stable time 


step far fourth order Range Kutta Method 


At-0.02min At=0.01 min 


T ime 
(min ) 

CPU 

X 

o 

CPZIxlO 7 

v ° 4 

10 

- 

- 

0.2175 

0. 1305 

20 

- 

- 

0.2193 

0.2607 

30 

- 

- 

0.2210 

0.3906 

40 

- 

- 

0 . 2228 

0.5201 

50 

- 

- 

0.2245 

0.6492 

60 

- 

- 

0.2263 

0.7780 


At=0. 

005 min 

At— 0.0001 min 

CPDxlO 7 

X a *1° 4 

CPUxlO 7 

X o x10 4 

0.83503 

0.1299 

0.83503 

0.1299 

0.8349 

0 . 26066 

0.8349 

0.2606 

0.8348 

0.39132 

0 . 8348 

0.3913 

0.8347 

0.5219 

0.8347 

0.5219 

0.8346 

0.65256 

0.8346 

0.6525 

0.83446 

0.7831 

0 . 83446 

0.7831 
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Iable 5 : Results of Isothermal Polymerization without gel effect 
(Results in paranthesis are numerical computations with At=0. 0001 min) 


Time [M]xl0 

-7 7 

CP ] xlO r 

Q 1 xl0 3 

°2 

X XlO 3 

O 


x 2 ,i<r 4 

PD I 

10 

0.8631 ' 

0.834 

0. 1745 

0. 7297 

0.1306 

0.3019 

0. 1386 

1.9851 


(0.8632) 

(0.834) 

(0. 1745) 

(0.7297) 

(0. 1305) 

(0. 3016) 

(0. 1385) 

(1.9864) 

20 

0.8340 

0.8328 

0. 1686 

0.6822 

0.2608 

0.5933 

0.2682 

1.9869 


(0.8341) 

(0.8328) 

(0. 1686) 

(0.6823) 

(0.2607) 

(O. 5930) 

(0.2680) 

(1.9873) 

30 

0.8059 1 

0.8317 

0. 1629 

0.6379 

0. 3907 

0. 8746 

0.3892 • 

, 1.9879 


(0.8060) 

(0. 8317) 

(O. 1629) 

(0.6380) 

(0. 3905) 

(O. 8741 ) 

(0.3890) 

(1.9882) 

40 

0.7788 

0.8305 

O. 1574 

0.5965 

0. 5203 

1. 1461 

0.5022 

1.9892 


(0.7788) 

(0. 8305) 

(0. 1574 J 

.(0.5966) 

(0.520) 

(1. 1454) 

(0.5019) 

(1.9895) 

50 

0.7526 

O'. 8294 

0. 1521 

0.5579 

0.6495 

1.4081 

0.6077 

1.9908 


(0.7527) 

(0. 8294) 

(O. 1521 ) 

(0. 5580) 

(0.6492) 

(1.4072) 

(0.6074) 

(1.9910) 

60 

0.7273 

0. 8283 

O. 1470 

0.5218 

0. 7784 

1.6610 

0.7063 

1.9929 


(0.7274) 

(0. 8283) 

(O. 1470) 

(0.5219) 

(0. 7780) 

(1.6599) 

(0.7059) 

(1.9930) 
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Table 6 : Results of Noni sotherma 1 Polymerization without gel effect 

(Results in paranthesis are numerical computations with At "0. 0001 m i n ) 

Real 

U = 2.0343x10 

hr cm c 


Dia = 2.0 cm 


Time Temp. CMlxlO 
(min) °K gmoles 


C P ] x 1 0 7 Q 1 xl0 3 
gmoles 


lit 1 it 


\ xlO 3 
o 



\ 2 xl0 4 PD 1 


10 331.32 0.8687 0.6631 0.1707 0.8791 

(331.49) (0.8662) (0.7529) (0. 1707) (0.7742) 

20 331.27 0.8462 0.6575 0.1662 0.8595 

(331.44) (0.8405) (0.749) (O. 1655) (0.7315) 

30331.23 0.8244 0.6531 0.1617 0.8013 

(331.39) (0.8159) (0.7456) (0. 1605) (0.6912) 

40 331.184 0.8034 0.6487 0.1575 0.7648 

(331. 342) (0.7920) (0.7423) (O. 1557) (0.6531 ) 
50 331.146 0.7830 0.6446 0.1534 0.7302 

(331. 296) <; 0.7609) (0.7392) (0.1510) (0.6174) 
60 331-110 0.7633 0.6407 0.1495 0.6772 

(331.252) (0.7467) (0.7361) (0. 1466) (0.5837) 


0. 09047 

0. 2454 

(0. 1142) 

(0. 2718) 

0. 1716 

0. 4714 

(0.2151) 

(0.5218) 

0.2515 

0. 6893 

(0.3191) 

(0.7749) 

0.3303 

0. 900 

(0.4221 ) 

(0. 1013) 

0.40809 

0. 1104 

(0.5242) 

(0. 1243) 

0.4849 

0. 1301 

(0.6254) 

(0. 1466) 


0.1329 1.9897 

(1.1332X1.9864) 
0.2575 1.982 

(0.2572) (1.9863) 
0.3755 1.9876 

(0.3738X1.987) 
0.4874 1.9874 

(0. 4835) ( 1. 9874) 
0.5936 1.9876 

(0.5868) (1.9883) 
0.6942 1.988 

(0.68403) (1.989) 



TABLE 7 


Effect of Shear on Rate of Polymerisation in Absence of 



Activated Carbon 


8.356 gmoles/lit 

P at 60° C 
m 

=. 0.869 gm/cc 

0.0533 gmoles/lit 

P at 60° C 

' nc 

= 1.0475 gm/cc 



R G" R R 


x 10 
1 i t/min 


4.28 


36.48 3.89 


5 

4 

Y x 10 

r p x 10 

0.5 

lit 

gmoles 

0.5 

gmoles min 

lit min 

8.466 

1.627 

i 

7.902 

1.514 



0.9995 


0.9985 


72.96 3.41 6.876 1.327 0.9917 


94.85 3.19 


6.407 1.245 0.9964 


120.39 3.55 


6.908 1.348 0.9988 


145.92 3.76 7.451 1.430 0.9935 


182.40 4.35 8.820 1.694 0.9991 













6 \ 




. u> r r 
r = 

r g~ r r 

sec" 1 

Slope 3 
x 10 

lit/min 

5 

Y x 10 

0.5 

lit 

0.5 

gmoles min 

O 

a 

o 

5.0 

10. 17 

36.48 

3.59 

7.267 

54.72 

4.19 

9.017 

72.96 

4.99 

10.37 




109.44 

4.63 

9.307 

131.33 

4.69 

9.565 


r p x 10 

gmoles 


lit min 



1.940 0.9947 


1.399 0.9995 


1.711 0.9997 


1.989 0.9850 


1.800 0.9992 


1.824 0.9968 










TABLE 9 


Effect of Activated Carbon (-350 mesh size) on 


Rate of Polymerization 


=8.356 gmoles/lit 
at 60 C = 0.869 gm/cc 


P at 60 C = 1.0475 gm/cc 
PS 


ClgU 



1.66 60 


Slope 

3 

x 10 

lit 

“< 

X 

mi, 

o 

* 

min 


3.868 

0.983 

4.582 

1.114 

4.300 

1.088 

4.560 

1.085 

4.520 

1.130 

5.220 

1.331 


x 10 P PO 

gmoles lit min 


lit min gmoles 
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TABLE 10. 

Effect of Activated Carbon (-250 mesh size ) on 
Rate of Polymerization 

CMD = 8.356 gmoles/1 it ^^at 60 * C = 1,0475 9 m/cc 

P at 60 C = 0.869 gm/cc 
m 


CCD 

o 

1 

ci 2 d 

Slope 

4 

Y x 10 

r P , 

1 

2 

R 

c 1 

x 10 

CCD 


3 

X 10 


4 

x 10 

r p “ r P0 


gmoles 

lit 

gmoles 

lit 


gmoles 

lit min 


lit 

gmoles 

lit 

min 


lit min 

gmoles 


0.00 

CO 

0.037 

3.868 

0.983 

1.581 

- 

0.997 

1.66 

60 

0.037 

3.960 

1 .009 

1 .618 

266666.7 

0.999 

2.00 

1 

50 

0.037 

3.952 

1.010 

1 .615 

289855.0 

0.993 

2.33 

42.86 

0.037 

3.822 

0.962 

1.562 


0.997 

3.00 

33.33 

0.039 

4.390 

1.087 

1.794 

46838.4 

0.999 

3.66 

27.27 

0.038 

4. 170 

1.047 

1.704 

80775.4 

0.999 
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* * ft * u f 


Gel and Glass effort * . „ ^ . 

~~ ct - const ltuti vp eq uations 


k t = 


k t« D 

to o 


D +k t e * *. 
o to t O 


k = k D 
P po o 


where 


D + k & \ 

0 Po p o 


(gel effect) 


(glass effect) 


& t ~ & ° e ><P(E dt /RT) 

CI 2 3 

% = & °p exp(E dp /RT) 

D o — r 2 . 303 ( 1 -^ ) 

[ L 

L A +B(i-y. ) J 
Y k tc = k ° ex P ( E /RT) 


(11.1) 

( 11 . 2 ) 

(11.3) 

(11.4) 

(11.5) 


( 11 . 6 ) 



